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Abstract 

We discuss several geometric PDEs and their relationship with Hydrody- 
namics and classical Electrodynamics. We start from the Euler equations of 
ideal incompressible fluids that, geometrically speaking, describe geodesies on 
groups of measure preserving maps with respect to the L metric. Then, we 
introduce a geometric approximation of the Euler equation, which involves the 
Monge- Ampere equation and the Monge-Kantorovich optimal transportation 
theory. This equation can be interpreted as a fully nonlinear correction of the 
Vlasov-Poisson system that describes the motion of electrons in a uniform neu- 
tralizing background through Coulomb interactions. Finally we briefly discuss 
an equation for generalized extremal surfaces in the 5 dimensional Minkowski 
space, related to the Born-Infeld equations, from which the Vlasov-Maxwell 
system of classical Electrodynamics can be formally derived. 
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1. The Euler equations of incompressible fluids 

The motion of an incompressible fluid moving in a compact domain D of the 
Euclidean space R"^ can be mathematically defined as a trajectory t —> g(t) on the 
set, subsequently denoted by G{D), of all diffeomorphisms of D with unit jacobian 
determinant. This space can be embedded in the set S{D) of all Borel maps h from 
D into itself, not necessarily one-to-one, such that 

(f){h{x))dx — I 4'{x)dx 

D Jd 

for all e C{D), where dx denotes the Lebesgue measure, normalized so that the 
measure of I? is 1. For the composition rule, G{D) is a group (the identity map / 
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being the unity of the group), meanwhile S{D) is a semi-group. Both G{D) and 
S{D) are naturally embedded in the Hilbert space H = L^(_D,R'^) of all square 
integrable mapping from D into R'^ and, therefore, inherit from H a formal Rie- 
mannian structure. The equations of geodesies on G{D) turn out to be exactly |AKj 
the equations of incompressible inviscid fluids introduced by Euler near 1750 jEu| . 
The Euler equations play a fundamental role in Fluid Mechanics (for geophysical 
flow modelling in particular) and their global well-posedness is one of the most chal- 
lenging problems in the field of nonlinear PDEs. Their mathematical importance 
is confirmed by the recent publication of several books by Arnold-Khesin )AK| . 
Chemin |Ch| . P.-L. Lions [Li|j Marchioro-Pulvirenti |MPj . as well as by Majda's 
lecture in the Kyoto ICM [CTa] , 

From a geometric point of view (different from the usual PDE setting which 
consists in solving the Euler equations with prescribed initial conditions), it is nat- 
ural to look for minimizing geodesies between the identity map and prescribed 
measure preserving maps. More precisely : 

Definition 1.1 Given h E G{D), find a curve t G [0, 1] — » g{t) e G{D) satisfying 
g{0) = I , g{l) = h, that minimizes 

Anig) ^\j\\g'{t)\\l.dt^\j'^j^\dtg{t,x)\'dxdt. 

The infimum is nothing but ■7^6^(1, h), where 6d denotes the geodesic distance 
on G{D), and any smooth minimizer g must be a smooth solution of the Euler 
equations (written in "Lagrangian coordinates" ) 

g" og-^ = -Vp, 

where p = p{t,x) G R is the pressure field and Vp = {dx^p, ...jOxaP)- The mini- 
mization problem will be subsequently called "Shortest Path Problem" (SPP). 
The basic local existence and uniqueness theorem for the SPP is due to Ebin and 
Marsden |EM| . If h and / are sufficiently close in a sufficiently high order Sobolev 
norm, then there is a unique shortest path. In the large, uniqueness can fail for the 
SPP. For example, in the case when D is the unit disk, h{z) = —z, the SPP has two 
solutions g{t, z) = ze"'"*'^* and g{t, z) = ze~*'^*, where complex notations are used. 
In 1985, A. Shnirelman | Sh| found, in the case D = [0, 1]^, a class of data for which 
the SPP cannot have a (classical) solution. These data are those of form 

h{xi,X2.,xz) = {H{xi,X2),x:i), 

where H is an area preserving mapping of the unit square, i.e. an element of 
G'([0,l]2), for which 

'5[o,i]3(^, ^) < '5[o,i]2(/,i?) < +00. 

(This means that, although h is really a two dimensional map, genuinely 3D mo- 
tions perform better to reach h from / than purely 2D motions.) 
Shnirelman also proved [SHI, EH], that 5([0, 1]'') is the right completion of G'([0, 1]'') 
for the geodesic distance 5, for all dimension d > 3. (Notice that S'([0, l]'^) is the 
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completion of G{[0, l]'^) for all d > 2 jNe| . So, the case d = 2 is very peculiar.) 
In such situations, a complete existence and uniqueness result for the SPP was ob- 
tained in |Br2| . provided the pressure field is considered as the right unknown and 
not the path t ^ g{t) itself. 

Theorem 1.2 Let h G 5([0,1]'^) of form h{xi,X2,X3) — {H{xi, X2), x^) with H S 
S'([0, 1]^). Then there is a unique vector-valued measure \/p{t, xi,X2) such that, 
for each sequence of curves t G [0,1] — > 5n(i) G ^([0,1]"^) labelled fey n e N and 
satisfying 

^[o,i]3(gTi) -> 2^[Qs]-'iI^h), ||g„(l) - ft.||i2([o,i]3 ^0, 
as n ^ 00, then (in the distributional sense) 

9n°9n^ -Vp. 

In other words, the acceleration field of all minimizing sequences converge to — Vp 
which uniquely depends on data h. The proof relies on an appropriate concept of 
generalized solutions (related to "Young's measures" |Yoj . |Taj . |DM| . |Shej ) that 
describe the oscillatory behaviour of the {gn) as n ^ +00 and reduces the SPP to 
a convex minimization problem. (See ,Br2. ^ for more details.) More precisely, the 
associated measures 

c„(t,x, a) = S{x -g„(i, a)), mn{t,x,a) = dtgn{t,a)5{x - gn{t,a)), 

have cluster points (c, m) that have the following properties : 

1) m is absolutely continuous with respect to c and its vector-valued density v{t, x, a) 
is c— square integrable; 

2) c and v do not depend on 2:3 and v^ = 0, 

3) c and v solve 

dtC+V,j:-{cv) = 0, 9f(cw) + Va;.(CT®w) +cVa;p= 0, (1.1) 

where the product cV^P has to be properly defined (in a way related to the work of 
Zheng and Majda [ZMp . Equations (|1.1() are obtained as the optimality equations 
of the convexified minimization problem. Therefore, it is a priori unclear they have 
any physical meaning as evolution equations. However, they correspond, up to a 
change of unknown, to the hydrostatic limit of the Euler equations, obtained from 
the Euler equations by neglecting the vertical acceleration term, namely: 

Kg"og-^^~Vp, 

where K is the singular diagonal matrix (1,1,0). These hydrostatic (or "shear 
flow") equations are widely used for atmosphere and ocean circulation modelling, 
as the building block of the so-called "primitive equations". However, they are 
more singular than the Euler equations and their mathematical analysis is very 
limited, as discussed in |^. Conditional well posedness and derivation from the 
Euler equations have been established in jBr3| and jGr| . 
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Remarks 

An intriguing question is whether or not the uniqueness of Vp can be proved 
by more classical tools even in the case when H £ G([0, 1]^) can be connected to 
the identity map by a classical shortest path on G([0, 1]^). 

Since S{[0, 1]^) is the right completion of G([0, 1]'^) with respect to the geodesic 
distance, one could expect the SPP to have a solution in S{[Q, 1]^) for all data h. 
This is not true. An example of such a data is h{xi, X2, x^) — (1 — xi, X2, x^). Only 
generalized flows, as discussed in |Br2) . can describe shortest paths in full generality. 

Example of generalized solutions 

Explicit examples of non trivial generalized shortest paths can be computed 
either numerically or exactly. Let us just quote a typical example, when D is the 
cylinder {{z,s) = (a;i,a:2,s), \z\ < 1, < s < 1} and h{z,s) = (— z,s). Then, 
the classical SPP has two distinct solutions g+{t, z, s) — (e"*z, s) and g-{t, z, s) = 
(e~*'^*2, s), with the same pressure field p = 7r^|zp/2, where complex notations are 
used on the disk \z\ < 1. (Notice that there is no motion along the vertical axis s.) 
Trivial generalized solutions are obtained by mixing these two solutions. However, a 
non trivial generalized solution exists and can be described as follows. For each fluid 
particle initially located at (z,s), the elevation s stays unchanged and the initial 
horizontal position z splits up along a circle of radius (1 — |zp)^/^ sin(7rt), with 
center 2Cos(7ri), that moves across the unit disk and shrinks down to the point —z 
as t = 1. In addition, each particle is accelerated by the pressure field p — 7r^|zp/2, 
as expected from the theory. 

2. Polar factorization of maps and the Monge- 
Ampere equation 

A way to define approximate geodesies on G = G{D) is to introduce a penalty 
parameter e > and to consider the formal (hamiltonian) dynamical system in the 
Hilbert space H = L'^{D,R'^) 

where the unknown M is a curve t M{t) e H, 6/5M denotes the gradient 
operator in ff, and 

d,(Af,G)= inf ||M-g|U (2.2) 

geG 

is the distance in H between M and G, where \ is the Hilbert norm of H. This 
approach is related to Ebin's slightly compressible flow theory |Ebj , and is a natural 
extension of the theory of constrained finite dimensional mechanical systems IRU| . 
jAK| . Notice that the approximate geodesic equation is sensitive only to the 
closure of G{D), which is, in the case D ~ [0, 1]"*, d > 2, the entire semi-group S{D) 
jNe) . As the penalty parameter e goes to zero, we expect that for appropriate initial 
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data, typically for M{t = 0) = / and {d/dt)M{t = 0) = uo, where vq is a smooth 
divergence free vector field on D tangent to the boundary, the time dependent map 
M converges to a geodesic curve on G. Because of the classical properties of the 
distance function in a Hilbert space, for each point M E H for which there exists a 
unique closest point tts{M) on S{D), we have 

Thus, we can formally write the approximate geodesic equation 1)2. l[l 

e^— M + M-^s(M) = 0. (2.4) 

Therefore, it is natural to address the following variational problem, that we call 
the Closest Point Problem (CPP) 

Definition 2.1 Given M e L^(D, R''), find h e S{D) that minimizes 

Jd 



2 



\M{x) - h{x)\^dx. 



The solution of the CPP is given by the Polar Factorization theorem for maps 



Theorem 2.2 Let M : D ^ R."^ be an L" map such that the probability measure 

Pm{x) — I 5{x — M{a))da 



JD 

is a Lebesgue integrable function on D. Then, there exists a unique closest point 
TTs{M) on S{D) and there is a Lipschitz convex function $ on TV^ such that 

7rs(M)(a) = (V$)(M(a)), a.e.aeD. 

In addition, $ is a weak solution, in a suitable sense, of the Monge- Ampere equation 

det{d^^<i>{x)) = pm{x). 

Thus, the Monge- Ampere equation |Ca| . which is usually considered as a non 
variational geometric PDE related to the concept of Gaussian curvature, also is the 
optimality equation of a variational problem closely linked to the Euler equations 
of incompressible inviscid fluids. In addition, the Polar Factorization theorem can 
be seen as a nonlinear version of the Helmholtz-Hodge decomposition theorem for 
vector fields which asserts that any L^ vector field on D can be written in a unique 
way as the (orthogonal) sum of the gradient of a scalar field and a divergence 
free field tangent to dD. Shortly after |Brl| . Caffarelli ,£a, established several 
regularity results for the Polar Factorization. For example, provided D is smooth 
and strictly convex, any smooth orientation preserving diffeomorphism M of D has 
a unique Polar Factorization with smooth factors, and tis{M) belongs to G{D). 
More recently, McCann |Mc| generalized the Polar Factorization theorem when D 
is a compact Riemannian manifolds. 
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3. Optimal Transportation Theory 

In fiii., , the solution of the CPP problem is based on the Optimal Transporta- 
tion Theory (OTT). The OTT was introduced by Monge in 1781 ^q) to solve an 
engineering problem and renewed by Kantorovich near 1940 |Ka| in the framework 
of Linear Programing and Probability Theory |RR| . In modern words, this amounts 
to look for a probability measure on a given product measure space AxB, with 
prescribed projections on A and B, that minimizes 



where the "cost function" c > is given on A x B. The CPP roughly corresponds 
to the case when A = B — D, c{x, y) — \M{x) — and each projection of ^ is the 
(normalized) Lebesgue measure on D. The connexion established in jBrl| between 
the OTT and the Monge- Ampere equation, enhanced by Caffarelli's regularity the- 
ory |Caj . introduced OTT as an active field of research in nonlinear PDEs. Let us 
first quote the work of Evans-Gangbo jEv| to solve the original Monge problem with 
PDE techniques, related to the Eikonal equations, and the recent contributions of 
Ambrosio, Caffarelli, Feldman, McCann, Trudinger, Wang. (A first attempt was 
made by Sudakov |Suj with purely probabilistic tools.) Let us next point out the 
importance of OTT for modelling purposes in Applied Mathematics. First of all, 
it is fair to say that the OTT and the Monge-Ampere equation were already key 
ingredients in CuUen and Purser's theory of semi-geostrophic atmospheric flows, 
which goes back to the early 1980s and preceded our Polar Factorization theorem 
(see references in |CNP| V Next, Jordan-Kinderlehrer-Otto , using OTT, es- 

tablished that the heat equation can be seen as a gradient flow for Boltzmann's 
entropy functional. More systematically. Otto Ot„ showed how the OTT confers 
a natural Riemannian structure to sets of Probability measures and recognized a 
large class of dissipative PDEs as gradient flows of various functionals for such 
Riemannian structures. Examples of such PDEs are porous media equations, lu- 
brication equations, granular flow equations, etc... Let us also mention that OTT 
has became a powerful tool in Calculus of Variations (through McCann's concept of 
displacement convexity |Mc p and Functional Analysis, where all kind of functional 
inequalities (Minkowski, Brascamp-Licb, Log Sobolev, Bacry-Emery, etc,...) can 
be established through OTT arguments, as shown, in particular, by Barthe |Ba| . 
McCann |Mcj . Otto, Villani "OV!. Let us finally mention that |BB| has provided 
for the OTT a formulation different from the Monge-Kantorovich one, by introduc- 
ing an interpolation variable (which was already present in McCann's concept of 
displacement convexity). This point of view is useful for both numerical 3B. and 
theoretical purposes, in particular, by allowing non trivial generalizations of the 
OTT related to section El 




AxB 



4. Approximate geodesies and Electrodynamics 
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Let us go back to the approximate geodesic equation (|2.4|) that can be (for- 
mally) written, thanks to the Polar Factorization theorem, 

duM{t,a) + {\/c^)it,M{t,a)) ^ 0, det{I - e^d^^c^{t,x)) = pM{t,x) (4.1) 

(where (l){t,x) stands for e~^(|a;|^/2 — ^{t,x))). A formal expansion about e = 
leads, as expected, to the Euler equation (written in Lagrangian coordinates) at the 
zero order and, at the next order (and exactly as d = 1), to 

duM{t,a) + iV^)it,Mit,a)) ^0, A^t, x) ^ 1 - puit, x), (4.2) 

which can be equivalently written as 

dtf + eV,/ - V^f^.V^/ = 0, e^A(l> =1- j fdi (4.3) 

by introducing the "phase density" 

f[t,x,C)^ [ S{x~ M{t,a))S{^~dtM{t,a))da. 
Jd 

This system is nothing but the Vlasov-Poisson system that describes the classical 
non-relativistic motion of a continuum of electrons around a homogeneous neutral- 
izing background of ions through Coulomb interactions. 

So, the approximate geodesic equation, which can be written as a "Vlasov- 
Monge-Ampere" (VMA) system, 

+ V,</).V^/ = 0, dct{I~e^d,,(b)^ J fdC (4.4) 

can be interpreted as a (fully nonlinear) correction of the Vlasov-Poisson system 
for small values of e. Recently, Loeper |Lo| has shown that the VMA system has 
local smooth solutions and global weak solutions. Loeper has also proved that the 
Euler equations and the Vlasov-Poisson system correctly describe the asymptotic 
behaviour of the VMA system as e ^ 0. The asymptotic analysis is based on 
the so-called modulated energy method already used in jBr5| to derive the Euler 
equations from the Vlasov-Poisson system. 

Notice that, thanks to the substitution of the Monge-Ampere equation (a 
fully non-linear elliptic PDE) for the classical Poisson equation, the "electric" field 
V0(i,x) is pointwise bounded by the diameter of D divided by e^, independently 
on the initial conditions. In particular, point charges do not create unbounded force 
fields as in classical Electrodynamics. 

5. A caricature of Coulomb interaction 

The approximate geodesic equation (|2.4(l can be easily discretized in space 
by substituting i) for D a discrete set of N "grid" points equally spaced in D, 
say Ai, Apf , ii) for H the euclidean space R'^^, iii) for G the discrete set of all 
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sequences {A^m ■■■,A„j^) G R'^-'^ generated by permutations a of the first A'' in- 
tegers, while keeping unchanged equation (|2.4|l . (Note that such a discretization 
using permutations cannot be so easily defined for the Euler equations, which for- 
mally correspond to the limit case e = 0.) Then M{t) = (Mi(t), ...,M]y{t)) can be 
interpreted as a set of N harmonic oscillators 

e'^M„ + M„-A,^(t) =0, (5.1) 

where the time dependent permutation a{t) is subject to minimize, at all time t, 
the total potential energy 

N 

J2\Mo.{t)-A,J\ (5.2) 

This system can be seen as a collection of N springs linking each particle Ma to one 
of the fixed particle Afs according to a dynamical pairing (3 = Ga it) maintaining the 
bulk potential energy at the lowest level. There is some ambiguity in the definition 
of this formal hamiltonian system for which the hamiltonian is given by 

Q— 1 Q— 1 

In particular, cT{t) is not uniquely defined at each time t for which several par- 
ticles have the same position. However, the potential is the sum of a quadratic 
and a Lipschitz concave functions of M . So its gradient has linear growth at in- 
finity and its second order partial derivatives are locally bounded measures. This 
is enough, according to recent results by Lions and Bouchut |Bo| . |Li2| . to ensure 
that unique global solutions are well defined for Lebesgue almost every initial data 
Ma(0), ^Mq(O), a = 1, ...,7V. As expected, the limit N +oo, e ^ (provided 
N goes fast enough to +oo), leads to the Euler equation, as proven in |Br4) . From 
the electrostatic point of view, the dynamical system describes a nonlinearly cutoff 
Coulomb interaction between N electrons (with positions Ma) and a background 
of N motionless ions (with fixed positions A^ ) ■ 

6. Generalized extremal surface equations and 
Electrodynamics 

As seen above, the approximate geodesic equation (|4.1|) — which has been in- 
troduced as a natural geometrical approximation to the Euler equations — turns out 
to be a model for electrostatic interaction with a non-linearly cutoff Coulomb poten- 
tial. This feature is somewhat reminiscent of the Born-Infeld non-linear theory of 
the electromagnetic field (Bij (see also B DLL| . |GZ| ...). Therefore, one may try to 
design from similar geometric ideas a non-linearly cutoff theory for classical Electro- 
dynamics. An attempt is made in jBr6| . Instead of considering springs linking two 
particles of opposite charges we rather consider (with a more space-time oriented 
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point of view) surfaces {t, s) X{t, s) spanning curves t X^{t) and t X^{t) 
followed by two particles of opposite charge, so that X{s = — l,i) — X-{t) and 
X{s = l,t) = X^{t), s G [—1,1] standing for the "interpolation" parameter be- 
tween the two trajectories. Just by prescribing {t,s) {t,s,X{t,s)) to be an 
extremal surface in the the 5 dimensional Minkowski space {t, s, xi, X3) (with 
signature ( — h + + +)), we get the building block of the model. In other words, the 
individual Action of each surface is 

J yTT^;^^^^^;TF"^^;T79;TFdtds, (6.1) 

(which is basically the Nambu-Goto Action of classical string theory). Next, we 
associate with X a "generalized surface" (p, J, E, B) (or more precisely a "cartesian 
current" in the sense of GMS ) defined by 

p{t,s,x)^5{x-X{t,s)), J{t,s,x)^dtX{t,s)5{x~ X{t,s)), (6.2) 

E{t, s, x) = dsX{t, s)5{x - X{t, s)), (6.3) 

B{t, s, x) = dsX{t, s) X dtX{t, s)5{x - X{t, s)) (6.4) 
and subject to compatibility conditions 

9,.p + V.S = 0, dtp + V.J^Q, dtE -dsJ ~y X B = 0. (6.5) 

In terms of (p, J, E, B) the Action of X can be written as 



Kip,J,E,B)= j ^p^-J^ + E^-B^. (6.6) 

Varying this Action under constraint (|6.5() leads to a system of evolution equations 
for {p, J, E, B) (see ^Br6, for an explicit form), that we can call "generalized ex- 
tremal surface equations" (GESE). They enjoy (at least in the simplest cases when 
the solutions depend on one or two space variables) many interesting properties : 
hyperbolicity, linear degeneracy of all fields BDLL , symmetries between t and s, 
J and E etc... From the GESE, we can derive through various (formal!) limiting 
process 1) the Born-Infeld and the Maxwell equations, as {p, J) are prescribed at 
s = — 1 and s = -|-1 (in which case there is no coupling between charged particles 
and the electromagnetic field), 2) the Vlasov-Born-Infeld and the Vlasov-Maxwell 
equations as {E,B) = is prescribed at s = — 1 and s = -fl (which corresponds 
to a free boundary condition dgX = for an individual surface and yields a full 
coupling between charged particles and the electromagnetic field). In spite of the 
possible physical irrelevance of the GESE, their mathematical analysis (global exis- 
tence, uniqueness, etc.), and the rigorous derivation from them of classical models, 
such as the Vlasov-Maxwell equations, are, in our opinion, challenging problems in 
the field of non-linear PDEs. 
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